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В.В. МКРТИЧЯН 


О РЕАЛИЗАЦИИ ПРОГРАММНОГО МОДУЛЯ СПИСОЧНОГО ДЕ- 
КОДЕРА СУДАНА ДЛЯ КОДОВ РИДА-СОЛОМОНА 


Решена задача разработки программного модуля списочного декодера Судана для 
кодов Рида-Соломона. Получена структурная схема и программная реализация де- 
кодера. Особенностью реализации является применение в модуле эффективного 
алгоритма факторизации Рота-Руккенштейна. 

Ключевые слова: коды Рида-Соломона, списочное декодирование, факторизация 
многочленов. 


Введение. Одним из мощных методов борьбы с помехами, искажающими 
информацию в системах ее передачи и хранения, является алгебраическое 
помехоустойчивое кодирование. Для применения этого метода необходимо 
выбрать подходящий код и эффективный кодек, работающий на этом коде [1]. 
В практике цифровой помехоустойчивой связи широкое распространение 
получили коды Рида-Соломона (РС-коды), при этом перспективным направ- 
лением остается разработка для них новых кодеков, способных эффектив- 
но работать при большей зашумленности канала. Большим достижением в 
теории помехоустойчивого кодирования было создание М.Суданом в 
1997 году принципиального списочного декодера для РС-кодов [2], кото- 
рый использует интерполяцию и факторизацию многочленов двух перемен- 
ных над расширением базового поля Галуа и способен с полиномиальной 
сложностью работать за пределами конструктивного кодового расстояния. 
Для решения алгоритмически трудной алгебраической задачи факториза- 
ции полиномов двух переменных над полем Галуа наилучшим средством в 
настоящий момент является алгоритм Рота-Руккенштейна [3]. Декодер Су- 
дана производит поиск всех кодовых слов, удаленных от пришедшего по 
каналу слова на расстояние Хемминга, не превышающее значения некоторо- 
го управляющего параметра / Следует отметить, что проблема разработки и 
технической реализации жестких и мягких списочных декодеров для РС-кодов 
и более общих алгебро-геометрических кодов актуальна и в настоящее время 
активно исследуется [4, 5]. 

Постановка задачи. Цель данной работы -— разработка и реализация про- 
граммного модуля списочного декодера Судана для компьютерной модели 
цифрового канала передачи данных, использующего в конструкции декоде- 
ра факторизационный алгоритм Рота-Руккенштейна. 
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Принципиальный декодер Судана и метод факторизации Рота- 
Руккенштейна. Пусть Е. — поле Галуа мощности а с фиксированным при- 
митивным элементом и/ Ех] — кольцо полиномов переменной х над полем 
Ра; Ра[х,У] — кольцо полиномов двух переменных хи унад полем Е. [6]. РС- 
код длины п=д и размерности К=а+1 можно определить как множество 
векторов {(р(х,), 2%), Р(»),..., Вх), где х = и/', Е [1], а р пробегает 
множество информационных полиномов, то есть полиномов из Еа[х] степе- 
ни не выше а. РС-код достигает границы Синглтона, поэтому кодовое рас- 
стояние равно л-4. Взвешенную степень монома ху определим как +4} а 
взвешенную степень 4ед.(О) полинома О из Е[х,У|], определим как 
наибольшую из взвешенных степеней входящих в него мономов. Далее 
мощность произвольного множества А будем обозначать через |А]. 

Алгоритм Судана включает два основных шага: шаг интерполяции, 
на котором по принятому из канала слову строится полином двух перемен- 
ных специального вида, и шаг факторизации, где данный полином разлага- 
ется на сомножители, по которым можно построить список. Он имеет поли- 
номиальную оценку сложности О(/?). 

Входными параметрами декодера Судана является некоторый 
управляющий параметр КЕ {0;1;...;п}) и параметры РС-кода: длина пл, раз- 
мерность кода К, примитивный элемент и/ поля Ро. Для корректной работы 
алгоритма необходимо выполнение ограничения: 

> 4|_/2@+ /а |- [4/2] 
(обоснование см. [2]). 

При декодировании на вход алгоритма подается пришедшее по ка- 
налу слово у=(у)-:" в виде сетки {(х, и)} -:", где х = и/!. Декодер Судана 
производит поиск всех кодовых слов в пределах Хемминговой сферы, цен- 
тром которой является у, радиусом — величина г = п- & Полученные слова 
называются связными и образуют искомый список. Величину г будем назы- 
вать радиусом локации списочного декодера. 

Выходом алгоритма является список всех информационных полино- 
мов Хх) (ЕЕ [х]), удовлетворяющих условию |{1 Км) = И}| > Е. Из [2] выте- 
кает, что этот список содержит истинное информационное сообщение. 

Алгоритм Судана: /* Вход: Е, 4, 6 и сетка {(ж, у),.. (ж, и) */ 

Шаг 0. Вычислить внутренние параметры г, /: 

т = [а/2] -1 и /=Г (2(и+1)/а )*? ] -1. 

Шаг 1. (Интерполяция). Найти любой полином ОЕ Е [х,У|, удовле- 
творяющий следующим условиям: 

1) деда (9) < т+1а; 

2) У! е [1], (м, И) = 0; 

3) О не равен нулю тождественно. 

Шаг 2. (Факторизация). Разложить О на неприводимые сомножители. 

Шаг 3. Выдать все полиномы Кх) такие, что (у - Кх)) является де- 
лителем О, причем Км) = и, по крайней мере в Ёзначениях Е [7] 

Вычисляемые на нулевом шаге параметры /и /7 гарантируют суще- 
ствование полинома двух переменных О, удовлетворяющего всем требова- 
ниям первого шага алгоритма Судана, в частности, обнуляющегося во всех 
точках входной сетки. 

Для реализации первого шага алгоритма Судан предложил искать 
полином © в виде 
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О(жУ=Х,-0 Ук-о""729 Чху. (1) 
Для нахождения коэффициентов 4 рассматривается однородная 
система линейных уравнений 
уно ко" (у) =0, КЕ {0;1;...;п}), (2) 
не имеющая в настоящий момент специальных методов решения, то есть 
решаемая, например, методом Гаусса за время О(/?). Разрешимость системы 
обеспечивают значения параметров /, 17, рассчитанных на нулевом шаге. 

Замечание. Для возможности получения списков разных объемов в 
ходе декодирования значение управляющего параметра ЕЁ можно менять с 
перезапуском последнего шага алгоритма. 

Алгоритм Рота-Руккенштейна рассматривает полином двух перемен- 
ных О(х,у) (Е Ех) [3] как полином переменной у с коэффициентами из 
Р[Х: О(Х,У) = Хо О®(Х) уи ищет в кольце Е [х] его полиномиальные кор- 
ни, называемые у-корнями. 

В ходе работы алгоритм строит дерево коэффициентов полиномов 
Кх(ЕЕ[Х)), для которых разность (у - Кх)) делит полином О(х,у). Приме- 
няя алгоритм Рота-Руккенштейна в списочном декодере Судана, мы полу- 
чаем все полиномы искомого списка (это вытекает из структуры полинома 
О, см. доказательство в [2]), а также другие сомножители вида (у - Кх)). 
Все ненужные сомножители отсеиваются путем проверки связности входя- 
щих в них полиномов ЖХ). 

Этот алгоритм реализован в виде рекурсивной процедуры, 
конструкция и корректность которой основана на следующей теореме. 

Теорема. Пусть Кх)=У.> ЁХ -— у-корень ненулевого полинома двух 
переменных О(х, у) над полем Е., А(х)=Х.= Ех", где Р0. ПустьМ(х,У)= х® Ох, 
У), Ч-\(жУ)=Мкх, ху+й) , где О(жу)=О(жу), 60, В — максимальное целое 
число, такое что х№ делит О(х, У). Тогда для любого Ь0 выполняется: 
О(х, В(>)) = 0, М(0, ^) = 0. 

Входные параметры алгоритма Рота-Руккенштейна: К определяет 
размерность линейного векторного пространства, в котором осуществляет- 
ся поиск у-корней полинома О; /-— степень у-корня, при которой на данном 
шаге алгоритма отыскивается коэффициент; массив $ — вектор коэффици- 
ентов у-корня. Вызывать процедуру необходимо с параметрами (О, К, 0). 

Алгоритм Рота-Руккенштейна: 
ргоседиге Весоп$гисЕ (О(х, у)=0 Е Е [х,У, КЕ 7 *, Е 2) 

/* Глобальный массив $[0,...,А-1] */ 

Шаг 1. Найти максимальное у (Е 2*.{0}), для которого О(х, У/х- 
полином двух переменных. 

Шаг 2. М(х, у):= О(х, У/х. 

Шаг 3. Найти все корни у (Е Е.) полинома одной переменной М(0, у). 

Шаг 4. Для каждого из ‘у выполнять{ 

Шаг 5. $[Л := у. 

Шаг 6. Если / = А-1 ‚то выдать $[0],....5[&1] 

иначе { 

Шаг 7. М`` (х, У := М(х у+У. 

Шаг 8. М`(х, у := М``(х ху. 

Шаг 9. Вызвать Весопгисе (М` (х, у), К, #1). 


16 


Раздел “Математика” 








Итак, выходом алгоритма Судана для РС-кодов с алгоритмом Рота- 

Руккенштейна в качестве второго шага служит список, содержащий в каче- 
стве своего элемента истинное сообщение. 
Структурная схема списочного декодера. На рисунке, приведенном 
ниже, изображена схема декодирования, принятого по каналу вектора у, и 
система входов кодека. На вход декодера подаются размерность Ки длина 
РС-кода п, зашумленный кодовый вектор у, радиус локации / помеченные 
на схеме как ВХ1, ВХ2, ВХЗ, ВХ4 соответственно. Декодирование происхо- 
дит в блоках Б1-Б7. Рассмотрим подробнее работу этих блоков. 


Декодер Судана Входы декодера: 
ав 

вх! — Вычисление т, / ое 
вх2-—ф| Б! у > вх: 
[т г в 













Построение матрицы М 


вх3— | 52 


Нахождениие 4»; 


МО =0 









По 


строение , 
ажУ-Вн Эно" нару 







Факторизация © в (РУ 
65 Оку) = 2 9“) 









Хх) 


Формирование списка х) 
Б6 


| список Дх) 


вх4. 


Структурная схема списочного декодера, схема входов 


Блок Б1 на вход получает К, л и вычисляет значения параметров 
декодера /и 7. 

Блок Б2 строит матрицу однородной системы (2), обозначенную на 
схеме буквой М. Построение 1-Й строки происходит путем вычисления коэф- 
фициентов (х)“(у) при неизвестной а ( (х, И)е<(х, у)}=") и присвоения 
их координатам строки. При этом порядок следования элементов строки 
несущественен, но в реализации проще использовать естественный поря- 
док, получаемый при раскрытии двойной суммы. 

Блок БЗ, решая матричную систему М© =0 методом Гаусса, где 
© - вектор-столбец коэффициентов 4 полинома (1), находит коэффициенты 
дю. 

Блок Б4 по полученным в БЗ коэффициентам дк строит полином (1). 
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Блок Б5 представляет реализацию рекурсивной процедуры факто- 
ризации полинома (1), основанную на алгоритме Рота-Руккенштейна. Дан- 
ная процедура в процессе работы строит дерево коэффициентов полино- 
мов, ветви которого образуют список, включающий искомый. 

Блок Б6 на вход получает радиус локации кодовых сообщений / ко- 
торый может служить, например, параметром обратной связи в том случае, 
если декодер является элементом схемы, последующие компоненты кото- 
рой могут потребовать сокращения объема списка, либо может принимать 
максимальное значение (см. выше). Б6 кодирует полиномы образованного 
на предыдущем шаге списка и, выбирая из полученных кодовых слов 
связные, вычисляет искомый список. 

О программной реализации. Рассмотрим некоторые аспекты программ- 
ной реализации структурной схемы, описанной выше. 

Арифметический процессор, служащий основой для вычислений, 
был реализован на базе динамической библиотеки теоретико-числовых ме- 
тодов М/ИпМТЕ-5_3_2 (см. например ПЁЫр://5Ноир.пе/пЧ/). Данная библио- 
тека позволяет производить вычисления в полях Галуа и содержит реали- 
зации алгебраических алгоритмов и структур, необходимых для работы ко- 
дека. Например, полиномиальный класс, который удобно использовать для 
представления истинных информационных полиномов списков в блоке Б6, 
класс матриц, который удобно использовать при построении матрицы М в 
блоке Б2, алгоритм Гаусса, используемый в Б3З. К сожалению, в данной реа- 
лизации библиотеки не предусмотрено представление полиномов двух 
переменных над полями Галуа, однако, наличие такого представления 
необходимо для реализации списочного декодера Судана. Для получения 
такой возможности в процессе реализации декодера был написан класс 
представления полиномов двух переменных на языке С++, использован- 
ный в блоках Б4, Б5 для представления полинома 0. 

Полином двух переменных О Е«[х,У] рассматривается как полином 
одной переменной ус коэффициентами из кольца полиномов Е«[»|, то есть 
как элемент кольца (Е [Х])[У|. Такое представление позволило реализовать 
полиномы двух переменных в виде динамического массива указателей на 
объекты классов полиномов одной переменной. Далее был построен на- 
следник этого класса, включающий реализацию процедуры факторизации 
полинома двух переменных Рота-Руккенштейна, необходимую для декоде- 
ра Судана. Благодаря специфике вышеописанного представления возмож- 
на эффективная реализация всех вспомогательных процедур класса факто- 
ризующегося полинома двух переменных, обращающихся напрямую к по- 
лям и методам класса-предка. 

Декодер реализован в виде С++ класса, включающего в себя шаг 
расчета параметров и шаг интерполяции алгоритма Судана. Перенесение 
шага факторизации в класс представления полиномов двух переменных 
над полями Галуа значительно снизило функциональную нагрузку на деко- 
дер. Эффективность реализации была также повышена путем классифика- 
ции управляющих параметров декодера на константные и варьируемые. В 
случае, когда при декодировании изменяются лишь варьируемый параметр 
г - радиус локации, - представляется возможным быстро получить искомый 
список на основе полученных ранее списков без перестроения структуры 
декодера. 

Структурная схема программно реализована под операционные си- 
стемы \Ипаом/5 95/98/МТ/2000/ХР на основе библиотеки классов МЕС [7]. 
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Выводы. Решена задача разработки программного модуля списочного де- 
кодера Судана для кодов Рида-Соломона с применением алгоритма факто- 
ризации Рота-Руккенештейна. При работе на низкоскоростных кодах дан- 
ный декодер способен исправлять большее количество ошибок по сравне- 
нию с другими декодерами [2, 4], но при этом вместо одного передаваемо- 
го сообщения выдает список сообщений, одним из элементов которого яв- 
ляется истинное. Поэтому для применения этого декодера в цифровых си- 
стемах передачи данных помимо разработанного модуля необходим блок, 
реализующий выделение истинного сообщения из списка. Проблема по- 
строения такого блока рассмотрена в работе [8]. Результаты настоящей ра- 
боты частично представлены в [7]. 
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Тре ргодгат тодше оЁ бидап д®егттеа 115 десодег Гог Веед-боотоп содез 
|5 аеу@оред. ТНе Ыоск Фадгат ап ргодгап геа|хаНоп оЁ Чесодег аге 
соп5гисеа. А рам оЁ {Ие геа|хаНоп 15 Базе оп {Ие еНесйуе ГасопгаНоп 
а!догйт ог Ко-Кискеп%еп. 
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